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We present an analytic universal impurity solver for strongly correlated electrons. We extend the 
many-body perturbation expansion via suitable two-particle renormalizations from the Fermi-liquid 
regime to the critical region of the metal-insulator transition. The reliability of the approximation 
in the strong-coupling limit is demonstrated by reproducing the Kondo scale in the single-impurity 
Anderson model. We disclose the origin of the Kondo resonance in terms of Feynman diagrams and 
find criteria for the existence of the proper Kondo asymptotic behavior in approximate theories. 
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I. INTRODUCTION 

One must have at hand a reliable theory capturing 
the salient features of the strong-coupling asymptotics to 
comprehend fully the effects of electron correlations in 

' metals. Most theoretical approaches are based on the 
well established weak-coupling expansion and we have 
only a few opportunities to test its extensions to the 
strong-coupling regime. The most reliable check of credi- 
bility of a theory in the strong-coupling regime is to com- 
pare its results with one of the existing exact solutions 
in this limit. There are two principal exact results for 

I strongly correlated electrons. It is the single-impurity 
Anderson model (SIAM) and the Id Hubbard model. 

[ The ground state of both models at half filling was con- 
structed with the aid of the Bethe ansatz. The former 
is a heavy- fermion liquid in the strong-coupling limit, ^ 
while the latter is an electron-hole liquid for arbitrary 
interaction strengthi^ 

Recently SIAM has won on importance, since its so- 
lutions form a fundamental building block in the con- 

[ struction of a dynamical mean-field theory (DMFT) of 
strong electron correlations via the limit to infinite- 
dimensional lattices.- Unfortunately, there are only ap- 
proximate impurity solvers in the strong-coupling regime. 
The only analytic theory, non-crossing approximation, is 
not a Fermi liquid and displays low-energy pathological 
features.^ Numerical quantum Monte-Carlo simulations'^ 

] are restricted to rather high temperatures. The most re- 
liable impurity solver proved to be the numerical renor- 

' malization group. It produces the expected three-peak 

1 spectral function with the central Kondo resonance for 
intermediate couplingi^ii Due to limitations of the nu- 
merical procedure it cannot be extended to arbitrarily 
large interactions. Neither of these solutions, however, is 
capable to produce a detailed ultimate picture of the the 
strong-coupling limit and to decide whether and in which 
form the Kondo resonance survives or whether the sys- 
tem undergoes a metal-insulator transition. This ques- 
tion can only be answered by an analytically controlled 
theory. 

There is presently no global analytic theory capable 
of tracing down the genesis of the correlation-induced 



Kondo resonance The standard weak-coupling, Fermi- 
liquid-based perturbation expansion renormalized with 
various self-energy insertions becomes unstable for in- 
termediate and strong electron interaction.^ Attempts 
to go beyond one-particle renormalizations by introduc- 
ing explicit vertex corrections and a two-particle self- 
consistency either do not lead to analytically solvable 
equations or have not yet produced the desired Kondo 
asymptotics in SIAM.i^iiS 

The aim of this paper is to construct an analytically 
controllable approximation that would reliably interpo- 
late between the weak- and strong-coupling regimes in 
models with a screened electron repulsion. To this pur- 
pose we use an appropriately renormalized diagrammatic 
many-body expansion. We demonstrate that it is not 
the one-particle self-consistency but rather a two-particle 
one that must be introduced in the critical region of the 
metal-insulator transition. We further show that only 
a symmetric inclusion of electron-electron and electron- 
hole scatterings leads to the proper critical behavior. 
The two-particle self-consistency will be achieved within 
the parquet approximation in which we separate singu- 
lar and regular functions. Not to lose analytic control 
only the potentially divergent functions with long-range 
fluctuations are kept dynamical while the regular ones 
with short-range fluctuations are replaced by constants. 
Such a simpliflcation does not affect the universal fea- 
tures of the critical behavior. We explicitly construct the 
approximation for SIAM so that its reliability becomes 
transparent by qualitatively correctly reproducing the 
exact weak-coupling (Fermi-liquid) and strong-coupling 
(Kondo) regimes. We identify with our construction the 
minimal necessary conditions on approximate schemes to 
cover the Kondo strong-coupling behavior. 



II. MODEL AND BASIC EQUATIONS 

We formulate the construction of an approximate the- 
ory interpolating between the weak and strong coupling 
regimes independently of the underlying model via renor- 
malized Feynman diagrams and equations of motion. To 
demonstrate reliability of our approach we, however, use 
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in our explicit calculations the single impurity Anderson 
model the Hamiltonian of which reads 

kcr 

We denoted rfi^ = dj-do-. When calculating the grand 
potential and thermodynamic properties of this impu- 
rity model we can explicitly integrate over the degrees 
of freedom of the delocalized electrons. To this pur- 
pose we standardly replace the local part of the prop- 
agator of the mobile electrons by a constant A(e) = 
TT^,^ |VkP<5(e — e(k)) = A the value of which we set as 
the energy unit. For simplicity we assume half- filled case 
IX = Ed + U/2. The impurity grand partition function can 
then be represented via a local Grassmann functional in- 
tegral 

Z = y Vi/jVip* exp 1^ ?/'n(*'^n + A + «sign(w„) A)7/;„ 
-U J^r n'^iT)nfir)^ . (2) 

For this integral we can straightforwardly build up the 
standard weak-coupling (diagrammatic) perturbation ex- 
pansion in powers of the interaction strength U (dia- 
grammatically represented as a vertex) and the bare 
propagator, being in the impurity case Gq{x + iy) = 
l/{x + i sign(?/)(A -I- \y\)) (diagrammatically represented 
as an oriented line). 

We know that in the particular case of SIAM the 
non-renormalized weak-coupling expansion is preferred 
and converges for arbitrary interaction strength. To en- 
able a treatment of phase transitions and singularities 
in more general extended lattice models we have, how- 
ever, to renormalize the bare expansion. The fundamen- 
tal quantity for renormalization in the many-body per- 
turbation expansion is the self-energy. In our approach 
we do not use the self-energy as the approximation- 
generating quantity determined directly from sums of se- 
lected (renormalized) diagrams. Instead, we represent it 
by means of the Schwinger-Dyson equation and the two- 
particle vertex F as follows 



So.(?a;„) = ^y^G'-g-(i^n') 



1 - E G<y{^^n+m) 



(3) 



where ujn = {'2n + 1)t:T and Vm — 2'KmT are fermionic 
and bosonic Matsubara frequencies in units fcs = 1, re- 
spectively. We assigned dynamical variables to the two- 
particle vertex as shown in Fig. [TJ 



FIG. 1: Diagrammatic assigning of dynamical variables, fre- 
quencies and spins, to the vertex r„.„.i ( 



We further introduce two-particle irreducible vertices, 
that is, vertices that cannot be diagrammatically dis- 
connected by cutting specific pairs of one-particle lines 
(propagators). The choice of the two-particle irreducibil- 
ity is ambiguous. Here we choose only the singlet 
electron-hole and electron-electron scattering channels. 
We denote the respective irreducible vertices A'^^ and 
A'^'^. With their aid we represent the full two-particle 
vertex via two nonequivalent Bethe-Salpeter (BS) equa- 
tions. The Bethe-Salpeter equation in the eh channel 
reads 



X Ti;i{iu!n" ,iuJn';ii^m) ■ (4a) 
In the ee channel we obtain 

- E ^f! ' ' «f^ni+n' )G't («W„" ) 

n" 

X Gl{iuJn+n'+m-n")'ri;l{iuJn" ,iuJn';i'^7n+n~n") ■ (4b) 

Equations ^ and ([4]) reduce the solution of the inves- 
tigated model to the determination of the irreducible ver- 
tices A'^'' and A'^'^. The principal idea of the parquet ap- 
proach is to use the topological nonequivalence of differ- 
ent Bethe-Salpeter representations of the full two-particle 
vertex via the irreducible ones. Due to the nonequiv- 
alence of different types of two-particle scatterings the 
reducible vertex in one channel is irreducible in the other 
oneJ^ If the vertex irreducible in both channels is the 
bare interaction U (parquet approximation) we obtain 
the basic parquet equation in our two-channel scheme 



Ti 



u 



(5) 



We use this representation in Eqs. ^ to exclude the full 
vertex F from the BS equations that then become cou- 
pled convolutive nonlinear integral equations for the irre- 
ducible vertices A'^'* and A^^ . These equations cannot be 
solved analytically and their numerical solution is avail- 
able only far from the Kondo regime. To gain an ana- 
lytic assessment of the low-temperature, strong-coupling 
behavior of the parquet equations ([U using Eq. ([5]) we 
have to resort to simplifications. 
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III. KONDO ASYMPTOTICS 
FUNCTIONS 



VERTEX 



right-hand side of Eq. ([5a|) . It reads 



A. Simplified parquet equations 

Equations (|4]) can be analytically solved (in an ap- 
proximate way) individually for the full vertex F by ne- 
glecting Eq. ([5|) and using the irreducible vertices as in- 
put. Approximations of this type with rcnormalized one- 
particle propagators are now called fluctuation-exchange 
(FLEX) . In the simplest case the irreducible vertices are 
replaced by the bare interaction. We find that vertex 
F tends to develop a pole in the BS equation with mul- 
tiple electron- hole scatterings, Eq. (Pa|) . caused by the 
imminent metal-insulator transition at zero temperature. 
While the two-particle vertex F remains finite (bounded) 
when calculated from Eq. (j4bp . We use the fact that 
only equation (|4ap develops a pole in simplifying the 
parquet approximation in the critical region of the metal- 
insulator transition to an analytically controllable theory. 
Since the output of the BS equation from one channel en- 
ters the kernel of the other one, only K'^'^ becomes singu- 
lar in the critical region of the metal-insulator transition. 
The irreducible vertex K"^^ remains bounded everywhere 
in the parquet approximation. 

We restrict our interest to the leading critical behavior 
of the solution of the parquet equations in the asymptotic 
region of the metal-insulator transition. To this purpose 
we use the reasoning of the renormalization group and ne- 
glect the short-range fluctuations. We take into account 
explicitly only long-range (in time) fluctuations shaping 
the singular behavior. In this way we do not influence 
the universal part of the critical asymptotics. We hence 
keep only singular functions dynamical (dependent on 
the relevant frequency only). Regular functions, where 
only short-range fluctuations contribute, can be replaced 
by constants, suitably chosen averaged values. We must 
further guarantee that only the leading low-frequency be- 
havior of singular functions determines the critical behav- 
ior in the Bethe-Salpeter equations. 

As a first step we approximate A'^'' with a static ef- 
fective interaction K^'^ = U. Inserting this ansatz into 
Eq. ([4a|) we obtain a FLEX-type equation for the irre- 
ducible vertex from the electron-electron channel 



1 + UXeh{Wm) 



(6a) 



where we denoted the electron- hole bubble Xeh{ii^m) = 
(3~^^^G^{iuJn)G i{iuJn+m)- Sincc the static approxi- 
mation to the electron-hole irreducible vertex K'^^ — U 
causes a deviation from the exact high-frequency limit 
^eh _^ jj ^ approximation ([5a|) holds in the low-frequency 
limit. We hence have to consider in our simplification 
only the leading low-frequency contribution from the 



A|j^ ^(il^m) = - = 



2 

U Xeh{iVrn) 



(6b) 



1 + UXeh{iVm) 

We now use Eqs. jibl and (O, and A^*"^ from Eq. ((6b)) 
to determine the effective interaction U . The right-hand 
side, of Eq. (|4bj) remains frequency-dependent even with 
our ansatz. Since we neglect all finite (short-time) fluctu- 
ations in regular functions, we have to replace the vertex 
F resulting from Eq. (j4bp with a constant. This replace- 
ment is not uniquely defined but the ambiguity has no im- 
pact on the qualitative (universal) critical behavior of the 
solution. We found that the most appropriate way in re- 
placing the full vertex with a constant is to multiply both 
sides of Eq. (|4b[) with pairs of one-electron propagators 
G{iuJn)G{iuJm-n) from the left and by G{iujn')G{iujm-n') 
from the right. We then average over the fermionic Mat- 
subara frequencies with a constraint m ~ 0, conserved 
during the electron-electron multiple scatterings. Doing 
so we obtain an explicit equation for the effective interac- 
tion, the static part of the electron-hole irreducible ver- 
tex, 



UXeeiO) = UXee - 



GfGl^^ 



(7) 



We denoted Xee — P ^ X^n (i'^n)Gx (iw-„) the static 
part of the electron-electron bubble and 



{G^G^X) = ^ E {iuj-n)X(iu^) . (8b) 



Equations (j6bp - ^ form a closed set of relations deter- 
mining the static effective interaction U and the dynam- 
ical vertex K""''^ {ivrn) as functionals of the one-particle 
propagators and the bare interaction U. Together 
with Eqs. ([5]) and ^ we have an analytic approxima- 
tion in closed form. The one-electron propagators may 
be either bare or renormalized with the self-energy, when 
one-particle self-consistency is used. The latter is not 
mandatory and we show that it worsens the approxima- 
tion with the Hartree one-electron propagators, in par- 
ticular in the intermediate and high-frequency sectors. 

The critical region of the metal-insulator transition 
in SIAM corresponds to the strong-coupling regime, 
C//A — + oo. It is reached when the denominator on 
the r.h.s. of Eq. (|6bp approaches zero. At zero tem- 
perature, half filling, and in the rotationally invariant 
case, G| = G| we define a dimensionless Kondo scale 
a = 1 + Uxeh{^) = 1-17 f_^dio^G+{ujf]/n ^ 0. 
We used an abbreviation G±{ijj) = G{uj ± iO+). The 
denominator of the r.h.s. of Eq. (|6b[) can then be ex- 
panded only to the lowest (linear) order in frequency if 
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we are interested in the leading singular behavior. We 
obtain 1 + Uxeh{t^ + lO"*") = a — iirUpQUJ. We denoted 
the density of one-particle states at the Fermi energy 
Po = — 7r~^3G'+(0). It does not depend at zero temper- 
ature and half filling on the interaction strength. This 
singular low-frequency asymptotics of A*™^ allows us 
to evaluate the leading contributions to the integrals in 
Eqs. ®. 

The leading contributions to the integrals with vertex 
\smg ^Yie asymptotic limit a — > read 



G{z)G{-z) 



In a\ 



I Inal 



|lna|\ f° duj 



f ^3[G+H2G_(-c.)2] , (9b) 



and 



(G^Gi^^) 



\lna\ 



duj 



3 [G+{ujfG-{-uj)^ 



(9c) 



The effective interaction has an asymptotic solution for 



UXe 



lna| ((GG)L) " 
^'pI {{GG)l,)_ 

{{GGr^J 



(10) 



where we denoted 

-^-i/"oorf^S[G'+M"G_(-a;)"]. 

The dimensionless Kondo scale measuring the dis- 
tance from the metal-insulator transition is determined 
from an equation for the critical point of the metal- 
insulator transition. It reads 1 + Uxeh{0) = 0. At 
half filling Xeh{0) = — Xee and we obtain a leading-order 
asymptotic solution for the distance to the critical point 



exp ■ 



2 2 \{GG)l^) 

Po ur<r<\s \ [^^e 



((GG)L) 



1] 



(11) 



Solution ([TTjl holds at zero temperature and half filling 
when the argument in the exponential tends to infinity 
and Uxee > 1- The latter condition can be viewed as 
a characteristic of the strong-coupling regime, since the 
FLEX, weak-coupling solutions are distinguished by the 
opposite inequality, Uxee < 1- 



B. Hartree propagators 

To manifest analytically that Eq. (fTTj) reproduces cor- 
rectly the Kondo asymptotics we use the bare (Hartree) 
one-particle propagators in Eq. pT|) . In this case the inte- 
grals with powers of the one-particle propagators can be 



exphcitly evaluated, {{GG)l 



r2(»- 



.l) (2n- 

'Po 



V(2n-1). 



Using this result we obtain an explicit representation for 
the Kondo scale 



a = exp<; [Upo - 1] 



(12) 



The exact result is a = exp{— 7r^t7po/8}. Determination 
of the exact prefactor at the linear dependence of the 
exponent of the Kondo scale on the interaction strength 
is beyond the reach of the present approximation. The 
prefactor depends on the way in which we replace reg- 
ular frequency-dependent functions with constants. The 
universal feature of the Kondo scale is linearity in the in- 
teraction strength of the exponent. This linearity is the 
major achievement of our construction. The other Fermi- 
liquid-based approximations fail in the strong-coupling 
regime of SIAM. The single-channel FLEX-type approx- 
imations either do not see the metal-insulator transi- 
tion at all (ee channel) or predict a quadratic exponent 
a ~ exp{— C/^Pq} {eh channels).— 



IV. KONDO ASYMPTOTICS - 
FUNCTIONS 



ONE-PARTICLE 



The core of the parquet approach is a self-consistency 
determining the two-particle irreducible vertices from the 
completely irreducible vertex (bare interaction in the par- 
quet approximation) and the one-particle propagators. 
Both, the completely irreducible vertex and the one- 
particle propagators serve as input to the parquet equa- 
tions. To complete the parquet approximation and make 
it conserving and thermodynamically consistent we have 
to determine the physical one-particle functions, that is, 
the physical one-particle propagator and the self-energy. 
All physical quantities in conserving theories are then 
generated from the full one-particle propagator or the 
self-energy and their dependence on appropriate external 
sources The physical one-particle propagator need not 
be identical with the " auxiliary" one used in the parquet 
equations for the two-particle vertices. 



A. Self-energy in the parquet approach 

The fundamental quantity for any thermodynami- 
cally consistent and conserving approximation is the self- 
energy. Knowing the full two-particle vertex from the 
parquet equations we determine the self-energy from the 
Schwinger-Dyson equation ([3]). It reduces in the criti- 
cal region of the metal-insulator transition to a simple 
FLEX-type expression 



^ f ^ U y-^ G^^{iuj„') 
^ ' Pifl + UXehi^'^n- 



') 



(13) 



We can analytically continue the sum over Matsubara 
frequencies to an integral with Fermi function. At zero 
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FIG. 2: (Color online) Density of states at zero tempera- 
ture, half filling, and various interaction strengths without 
one-particle self-consistency. 

temperature the real part of the self-energy reads 
3fiE+(w) 



= -U 



e{uj) / + e{-uj) 

■oo 

° dx 



dx ^ G+{x + uj) 



TT 1 + Ux+ix) 



f^'^ dx 1 

-Ue{-u) -^G+{x + u;)^ 

Jo 1 + Ux+ix) 

The imaginary part has a representation 



(14a) 



9E+(tj) = -U 



LU) / + 0{-uj) 



dx 

TT 



3 G+{x + uj)^ = 



1 



l + Ux+ix) 



(14b) 



We distinguished with a subscript -t- the way the real 
frequency axis is approached in functions with complex 
variables. That is, G+{x) = G{x + iO+), x+{x) = 
Xeh{x + iO+). The onc-particle propagators on the r.h.s. 
of Eqs. (fT4|) may either be bare or full propagators with 
the self-energy from the left-hand side. Both construc- 
tions lead to conserving approximations when all thermo- 
dynamic quantities are determined from the self-energy. 

We first use the bare (Hartree) one-particle propaga- 
tors in Eqs. (fT4|) . The density of states calculated from 
the one-particle propagator containing the self-energy 
from Eq. ((M]) is plotted for various interaction strengths 
in Fig. [21 The central Kondo quasiparticlc is well formed 
and its width shrinks with increasing interaction strength 
in accord with Eq. See Fig. [3] for a detailed de- 

pendence of the quasiparticlc peak on the interaction 
strength. Not only the central quasiparticle peak but 
also the satellite Hubbard bands are well formed. Even 
though our approximation was justified for small frequen- 



FIG. 3: (Color online) Detailed dependence of the quasipar- 
ticle peak in the density of states on the interaction strength. 



cies, the parquet equations with the Hartree propaga- 
tors reproduce the high-frequency features of the spec- 
tral function surprisingly well. The Hubbard bands are 
positioned very close to the exact values of the atomic 
levels of the impurity electrons at ±[//2. 



B. One-particle self-consistency 

Usage of the Hartree propagators may seem inferior 
to a theory with renormalized one-particle propagators in 
the parquet equations. In fact, the opposite is true, even 
though the one-particle self-consistency adds new self- 
energy diagrams to those explicitly taken into account 
in the parquet construction of the self-energy from the 
two-particle vertex. The solution with renormalized one- 
particle propagators loses most of attractive features ob- 
tained in the parquet approach with the Hartree propaga- 
tors. First, the asymptotic formula for the Kondo scale, 
Eq. (jlip cannot be evaluated explicitly when renormal- 
ized one-particle propagators are used. Hence, the exact 
strong-coupling asymptotics cannot be reached analyti- 
cally. Second, the satellite peaks are completely washed 
out and merge with the central quasiparticle peak of the 
density of states as shown in Fig. [4l The central peak 
broadens with respect to the solution with the Hartree 
propagators. The solution with renormalized one-particle 
propagators shows overall worse agreement with the ex- 
act Kondo behavior derived from the Bethe ansatz. An 
analogous behavior was observed in a static simplification 
of the parquet equations. The failure of one-particle 
self-consistent theories to reproduce reliably the Kondo 
strong-coupling regime in impurity models is in accord 
with early analyses of SIAM using perturbation expan- 
sion in the interaction strength. 

The density of states resulting from the parquet 
approximation with renormalized one-particle propaga- 
tors resembles the solution of the renormalized RPA of 




FIG. 4: (Color online) Density of states for one-particle non- 
self-consistent (NSC) and self-consistent (SC) solutions com- 
pared. 



FIG. 5: (Color online) Real part of the self-energy for one- 
particle non-self-consistent (NSC) and self-consistent (SC) so- 
lutions compared. 



Suhl.l" The central quasiparticle peak is, however, much 
broader than in the FLEX-type approximations. It is 
even broader than that from the parquet solution with 
the Hartree propagators. This indicates that the criti- 
cal region of the metal-insulator transition is reached, if 
ever, for much stronger interactions in the solution with 
one-particle self-consistency than without it. The high- 
frequency behavior of the parquet approximation with 
renormalized one-particle propagators and the FLEX so- 
lution seem to behave similarly. One can prove analyt- 
ically that the one-particle self-consistent theories with 
simplified one- frequency two-particle vertex functions be- 
have in the high-frequency region universally and no 
Hubbard satellite bands emerge. 

To assess the high-frequency behavior of one-particle 
propagators we evaluate explicitly the leading-order con- 
tribution to the self-energy in the critical region of the 
metal-insulator transition, that is, a 0. We obtain 
from Eqs. (fT4)) 



5Ri;+(w) 



u 



arctan 



In a I 5RG+(w) 



In 



a 

—2 



3G+(tj) 



(15a) 



1 



(15b) 



If we now assume uj 3> aA we can further simplify the 
expression for the self-energy to 



S+(c.) 



C/|lna| 



n I \ , 7rsign(w) ~ ' 
2 ma 



(16) 



We decorated the high-frequency one-particle quantities 
with tilde to distinguish them from the full ones. For 
tilde functions we can introduce a dimensionless variable 



;7|lna| 



(17a) 



We resolve the one-electron tilde propagator by using the 
Dyson equation of SIAM. We obtain an explicit solution 



U\\na\ 



-Ml 



T 



(17b) 



We see that one-particle self-consistent theories lead to 
a universal high-frequency behavior with a semi-elliptic 
density of states spread over a large interval of order 



U\ lna|/C/7r2p2. The FLEX solution differs only in that 

U — U. The semi-elliptic form of the density of states in 
high frequencies is hence universal for one-particle self- 
consistent theories. The Hubbard satellite bands are 
completely smeared out. The reason for the nonexis- 
tence of the satellite peaks lies in the behavior of the 
real part of the self-energy in the low-frequency region. 
We compared the non-self-consistent and self-consistent 
solutions in Fig. [5] We can see that the sharp low- 
frequency structure of the real part of the self-energy is 
significantly smoothed and the peaks are broadened by 
the one-particle self-consistency. Most importantly, the 
height of the peaks near the Fermi energy is so much low- 
ered that |a;| > |3fiE(a;)|. Only if this condition is broken, 
as is the case in the non-self-consistent solution, satellite 
peaks emerge. 

V. CONCLUSIONS 

We presented a construction of a universal analytic 
impurity solver reliably interpolating between the weak 
(Fermi liquid) and strong (Kondo) coupling metallic 



7 



regimes. The approximation is a simplification of two- 
channel parquet equations where the bounded irreducible 
vertex is replaced by a static effective interaction and 
only the low-frequency singular vertex is kept dynam- 
ical. Such a procedure is justified in the critical region 
of the metal-insulator transition where the low-frequency 
asymptotics of the singular two-particle vertex becomes 
dominant. We set with our construction the minimal 
necessary conditions that must be fulfilled to reproduce 
the Kondo exponential scale. The Kondo behavior is ob- 
served in an approximate solution of SIAM if the low- 
frequency singularity in the two-particle vertex domi- 
nates the dynamics and when electron-hole and electron- 
electron multiple scatterings are self-consistently mixed 
in a balanced way. 

We analyzed two versions of the parquet approxima- 
tion. In the first one, we used the Hartree one-electron 
propagators in the parquet equations for the two-particle 
irreducible vertices and the Schwinger-Dyson equation 
for the self-energy. In the second one we used the fully 
renormalized one-particle propagators. Although the lat- 
ter construction contains more Feynman diagrams renor- 
malizing the self-energy, its results are less reliable than 
the results of the former approach. This conclusion is 
not new and is also understandableJ^ii^ Renormaliza- 
tions due to the one-particle self-consistency unrealisti- 
cally smear and unfold the low-frequency structure of the 
self-energy. Consequently, the high-frequency features of 
the one-particle propagator are washed out and no Hub- 
bard satellite bands can be observed. 

One-particle self-consistency is standardly demanded 
in order to guarantee conservation laws and thermody- 
namic consistency. It is, however, not a necessary con- 
dition for approximations to be conserving. Conserva- 
tion laws are guaranteed in any theory with an approx- 
imate self-energy functional containing physical external 
sources. These sources are then used to generate the de- 
sired physical quantities via linear response theory. On 
the other hand, theories aiming at a description of crit- 
ical phenomena such as magnetic phase transitions or 
a metal-insulator transition must be self-consistent in 
some way. Self-consistency is necessary for any theory to 



handle singularities. Since singularities in models with 
itinerant electrons emerge only at the two-particle level 
in Bethe-Salpeter equations, we must introduce a two- 
particle self-consistency. The parquet construction of- 
fers a very natural way to reach this goal and to treat 
singularities in the Bethe-Salpeter equations appropri- 
ately. One-particle self-consistency influences the critical 
behavior of the two-particle vertex, but it does not offer 
any direct control of singularities in the Bethe-Salpeter 
equations. It can neither guarantee integrability of singu- 
lar vertices nor is capable to reproduce the correct critical 
behavior in the Kondo regime. When used in the par- 
quet approach, the one-particle self-consistency interferes 
in the control of the two-particle singularities from the 
Bethe-Salpeter equations achieved by the two-particle 
self-consistency. The one-particle self-consistency can 
hence be used only when appropriately compensated so 
that it does not significantly affect the two-particle crit- 
icality. This is, however, not the case in the parquet 
approximation and the non-self-consistent one-particle 
propagators deliver better results than the self-consistent 
ones. 

To conclude, we derived a global approximation that 
is analytically tractable, sufRciently simple, and univer- 
sal. It can be used in a number of physically interest- 
ing situations, including realistic (multi orbital), material 
specific models as kind of a mean-field theory for models 
with strong local electron interaction. It may stand as an 
alternative to a recently proposed impurity solver for the 
strong-coupling regime, Unlike the solver from Ref. [l^ 
our approximate scheme reproduces correctly the Kondo 
strong-coupling regime. 
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